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Abstract. In this paper we analyze the spin- wave excitations (magnons) of an 
inhomogeneous spin system within the Boltzmann-Gibbs framework and then connect 
the results with the nonextensive approach (in the sense of Tsallis statistics). 
' Considering an equivalence between those two frameworks, we could connect the 

entropic parameter q with moments of the distribution of exchange integrals of the 
^ inhomogenous system. It ratifies the idea that the entropic parameter is connected to 

jH . the microscopic properties of the system 
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1. Introduction 



o 

(N 

^ ! Inspired in multifractals, Tsallis proposed a generalization of Boltzmann-Gibbs entropy 

: (Sbg) [i\ 

q. S, = k (qe^) (1) 



where q is the entropic parameter for a specific system and is connected to its dynamics 



I as recently proposed [21 [3], are the probabilities satisfying J2i Pi = ^i k is a constant, 
and limg^i Sg = Sbg- This entropy for a system composed of two independent parts A 
and B, such as the probability is given by P{Af]B) = P{A) P{B), has the interesting 
property of nonextensivity (see for example Refs. ^E]) 

Sg{A + B) = Sg{A) + Sg{B) + (1 - q)Sg{A)Sg{B) (2) 

Besides representing a generalization, Sg, like Sbg, is nonnegative, concave, Lesche- 
stable (Vg > 0). Recently it has been shown that it is also extensive for some sorts of 
correlated systems [6]. 

Tsallis statistics, or nonextensive statistics, attempt to handle some anomalies 
that appear in physical problems which cannot be treated with Boltzmann-Gibbs 
(BG) statistics, for instance, long-range correlations, intrinsic cooperativity, multifractal 
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structure, dissipation in mesoscopic scale, strong non-Markovian microscopic memory, 
etc [7]. These anomalies have the common characteristic of presenting power-laws, 
instead of the ordinary exponential-laws, which is also a characteristic of some complex 
systems. Its applicability ranges from protein folding [8] to financial markets [9], and 
from turbulence [10] to cosmic rays [H]. For example, in condensed matter problems 
we can cite Ising ferromagnets. Landau diamagnetism, electron-phonon systems, tight- 
binding-like Hamiltonians, metallic and superconductor systems [12]. In adittion, an 
interesting example is in [131 El where the authors could predict some peculiar magnetic 
properties of manganites using nonextensive statistics like nanoscale inhomogeneity and 
phase coexistence, fractal structures, and long-range interactions [IHl [16], [17] . 

Herein, we present some results comparing an inhomogeneous spin system within 
BG framework and a homogeneous spin system into a nonextensive approach. This 
comparison lead us to a connection between the nonextensive parameter q with specific 
moments of the distribution of the exchange integral of the inhomogeneous system. 
Thus, the spin waves in a inhomogeneous magnetic media can be described by the 
nonextensive statistics and the entropic parameter is connected to the microscopic 
properties of the system, as previously shown for other systems by Beck [18], Beck 
and Cohen [19], Wilk and Wlodarczyk [20], Reis et al. [2], and therefore can be seen as 
a measurement of its complexity. 

2. Spin- waves 

2.1. Magnons within inhomogeneous medium: Boltzmann-Gibbs framework 

In a ferromagnet at T = OK all the spins have the maximum projection S along the z 
direction; this is the ground state configuration [211 ES] . Letting the spin system be in 
thermal contact with a reservoir, as the temperature increases, it will leave its ground 
state, the projections along the quantization direction will be reduced, and a wave- 
like perturbation will flow through the spin system; that is the spin- wave (magnons). 
The spin-wave theory leads to the description of the magnetism of ferromagnets at 
low temperatures, in the regime where the total angular momentum is close to its the 
projection along the z direction, S = S^. 

We will consider a system of spins, each one interacting with z neighbors in a 
inhomogeneous way, and in the presence of a magnetic field Bq. Thus, the Hamiltonian 
for this inhomogeneous magnetic system is given by [211 1231 [21] 



in which Jd(R) > (always ferromagnetic) describes the inhomogeneity of the media, 
i.e., there is a distribution of exchange interactions /(J), and h = g hb Bq plays the role 
of the applied magnetic field. The whole Hamiltonian can be re-written in terms of the 
collective motion operators in order to give us the magnetization per unit of volume. 




(3) 



R,d 



R 
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One can write the spin operators as 

Sr • Sa+d = 2 '^'R.'^'R.+d + ^n^ii+d + S^^zS^+dz- (4) 

The Holstein-Primakoff transformation of spin operators, at low temperatures, is 
given by 

^ V2San and ^ V2S (5) 



in which and obey the commutation relation or, 
and Or can be written in terms of the collective motion of the system 



an 



J^e'^'^bk and 



1. The operators or 
(6) 



k V J'* A; 

in which [bk, b^] = 1. 

Thus, one can rewrite the exchange term of the Hamiltonian (j3]) 

Jd(R) SR-SR+d = -- ^ Jd(R) 'S'RS'R^^d + 'S'R'S'R+d] - X! ^rf(^) '^R-2'^(R-+d)z-C^) 

R,d R,d R,d 

Using the low-temperature Holstein-Primakoff transformation, in terms of the 
collective motion operators, we have 



1 



2S 



o E U^) ^R^R+d + SnS^^d = E I ^ E U^) cos(k • d) I 6+ b. 



" R,d 

and 



R,d 



R,d 



R,d 



k- 



(8) 



(9) 



R,d 



Note that we have excluded the magnon-magnon interaction, represented by the term 
of Uk rik. 

In terms of these operators, 

Srz = S - a^ttR (10) 

where S is the spin value per site and therefore the second term of the Hamiltonian can 
be written as 



hY.Siiz = hY,bt bk-hNS 



R 



where is the number of sites. Hence, the Hamiltonian becomes 

l-L = -{hNS + NS^J) + Y. (/^ + 25J-^E Jd(R)cos(k-d) j n^. (12) 

k \ ^ R,d / 

where = b'^bk is the Boson number operator and J = -^Y.-R.,dJd{^- The first 
two terms represent the fundamental state of the system, or the total energy without 
excitations. The term that describes the magnons is the second one. It has the form 
J^k^^k^k and gives the dispersion relation for this inhomogeneous magnetic system 



nujk = h + 2S J -"^Y. JdC^) cos(k-d). 

^ R,d 



(13) 
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For large wave length, one may write 

(14) 



^ Jd(R) cos(k-d)^^ Jrf(R 

R,d R,d 



l--(k-d)2 



and (k ■ d)^ = A;^ a^, being a the lattice parameter. Thus, the dispersion relation is 

huk^ h + k^V{J). (15) 

where T>{J) = a"^ z S J is the stiffness parameter and z is the number of first neighbors. 

Since the interaction varies between spins, one may consider that it has a 
distribution /(J). Thus, the average magnetization varies with respect to the saturation 
value of the magnetization and its variation is given by 

(Am) = 1^ rdJf{J) rdkk'{n,)j (16) 
where {nk)j is the Planck distribution. Then, flTBl) becomes 



4 7r2 Jo ' ' 'V{Jf/^ Jh/ksT 

where x = {k'^V{J) + h)/kBT. 

For i?o = 0, i.e., h = 0, the inner integral becomes 

X"*"/^ da; ^/7^ 



1 V^(3/2) (19) 
e^ — 1 2 

where ({n) is the Riemann Zeta function. The volume magnetization variation due to 
the magnon excitation of a inhomogeneous is then given by 

- ^ mTi:"^^ - ^ mr <-='^>- 

It is important to emphasize that the volume magnetization change of the 
inhomogeneous system has a T^/^ dependence (like the homogenous case) and also 
depends on the —3/2 moment of the distribution of exchange integrals (J^'^^^). This 
exponent is expected since 3 is related to the dimension of the system and 2 refer to the 
dynamics, i.e., came from the dispersion relation (fTSl) . 



2.2. Magnons within homogeneous medium: Nonextensive framework 

The dynamics of a system is given by its Hamiltonian Ti and the wave number fc, defined 
by 7i, is, consequently, related to the dynamics. On the other hand, the statistics of a 
system is given by an average over a great number of variables; and lies, for instance, 
in the number of Bosons for each wave number k. This average over weighted 
states makes it possible to obtain the relation of microscopic physical properties and 
macroscopic thermodynamic quantities such as the volume magnetization variation. 
The nonextensive approach proposes a change on the statistics of the system, not on 
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its dynamics. Thus, we assume an equivalent Hamiltonian as ([3]), but homogeneous in 
this framework, i.e., the exchange integral can be taken out of the sum. The dispersion 
relation is therefore given by = S J k'^ , where J is the exchange integral of this 
homogeneous system. The volume magnetization in this nonextensive scenario can be 
written as [25] 



(Am), 



2tx 



dk k ( rik) 



q,J 



(21) 



in which Uk = bk is the Boson number operator and (...)g is not the standard Planck 
distribution, but its g- version, i.e., generalized Plank distribution 

Tr{p«} E~,.„|i-(i-<;)(/3n...)l* ■ ' ' 

Using the dispersion relation described above and making (pT]) dimensionless, one 



gets 
{Am)g 
where 
f{x,q) - 



gfJ'B ( ksT 
47r2 \a^SJ 



3/2 



dxx^^'^ f{x, q) 



J:Z=o^k [l-(l-g)(xnfc)]T^ 
[l-(l-q)(xnk)]^ 



Finally, the magnetization can be written as 



(Am), 



9f^B I ksT 
47r2 \a^SJ 



3/2 



F{q) 



(23) 



(24) 



(25) 



where F{q) is the integral which appears in Eq. (!23l) . One can see that the magnetization 
in this scenario has the same T^/^ behavior as in (I2U]) . It is a consequence that neither the 
dynamics (e^ oc k"^) nor the dimension {d = 3) have being changed. All the information 
about the homogeneity and/or inhomogeneity of the system is into the statistics and, 
consequently, into the coefficient of the magnetization change. 

An analytical connection between the entropic parameter q and the volume 
magnetization variation can be obtained at the limit (g — 1) — > 0. At this limit, we 
can write (123]) as 



J 



Thus the volume magnetization change is now given by 



(Am), 



9J^ 
2 7r2 



dk k'^{nk) 



9[^B 



knT 



3/2 



47r2 \a^zSJ 



{q 



in which F' and F" are dimensionless integrals 



X 



^/^dx 



TT 



el' 



1] 



1 _L ^2 [£""+11 



5.2277 (g - r 



(26) 



(27) 



f28) 
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and 



r 



/„ |J-1|3 h^fclll. = 16.4154 -68.6515 to- 1). (29) 



As it is an approximation for q close to 1, there is no necessity of terms higher than 
' 1), so the second term of V"^ can be neglected. Thus, 



(Am), 



(30) 



47r2 {a^zSjJ 

One can see that, in this approximation, the magnetization is directly related to the 
entropic index and when g — 1, it recovers the usual result (homogeneous case within 
BG statistics). 

3. Mean field approximation and the critical temperature 



Let us consider that the two systems discussed before have the same critical temperature, 
as already done in the literature [Hj. Considering the Hamiltonian ([3]) for an 
inhomogeneous magnetic system within the mean field approximation, one may change 
the quantum operator Sa+d for its thermal average (SR+d)T- Thus, considering z first 
neighbors of an atom on the R^^ site of the lattice, its Hamiltonian becomes 

7^f, = -(j)^SR-(SR+d)T (31) 

where in this approximation we can consider the exchange interaction between the spins 
as an average value (J). It is reasonable, because all of the spins, in the mean field 
approximation, interacts with all others spins in the same way. 

For the above Hamiltonian it is straightforward to obtain the critical temperature 

(32, 

An analogous calculation can be done in the nonextensive scenario [26]. The 
generalized Brillouin function [2] gives us the critical temperature 

(33) 

3kB 

in which J' is the exchange integral in this framework. The relation between these two 
temperatures is given by [Ml [26] 

rj«) = grw. (34) 

Thus, using ( l32l) . ( l33i) and ( l34l) one find the relation between the two exchange 
integrals 

J={J). (35) 

that is, the exchange integral in nonextensive framework is equivalent to and average of 
the inhomogeneous one. This result is expected since, as already discussed above, we are 
not changing the dynamics of the system, but only but only the statistical treatment, 
which is used to calculate its thermodynamical properties. 
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4. Equivalence between the two frameworks 

Comparing tlie magnetization change per unit of volume in the inhomogeneous 
framework fl2T|) with its analogous in the nonextensive scenario (l25l) . one finds that 

^^^^ - ^H2j -~^\2) Jjy^ ^^^^ 

The above equation is a relation of the q parameter and moments of the exchange 
interaction distribution /(J). Figure ([1]) presents the expression above numerically 
solved for q G [0.1, 1.9]. This procedure comparing the magnetization was already been 
done p^, where those authors have found similar results inspired in Superstatistics |19j . 




0,5 1,0 1,5 2,0 

Entropic parameter q 



Figure 1. The entropic parameter q is connected to specific moments of the 
distribution of exchange integrals. This result is valid for any /(J) and shows that the 
entropic parameter is connected to the physical properties of the system. 
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An analytical connection between the entropic parameter q and the specific 
moments of the exchange integral of the inhomogeneous magnetic media can be obtained 
using the expression for the volume magnetization change in the limit (g — 1) — >■ 0. 
Comparing the (!2T|) and (!28|) one gets 



2 ^ 

or 



2 ) + 2.98 (g 



2 ^ 



(J)-3/2 



0.78 



(37) 



(3^ 



(J) -3/2 

The result above is also valid for any /(J) and shows that the entropic parameter 
is connected to the physical properties of the system [2|, [181 [13 120] . 



5. Final remarks 



Summarizing, in the present work we have shown that the q parameter can be seen as a 
measurement of the inhomogeneity of a magnetic system. It ratifies previous works [21 [3] 
in which, inspired in Superstatistics [19], the entropic parameter q was related to the 
first and second moments of the distribution of magnetic moments of manganites 

q{2-qf=^-^^ (39) 

and was also experimentally verified. Thus, the present work corroborates the idea 
that changing the usual Boltzmann— Gibbs statistics to one that is able to describe 
power-laws (Tsallis statistics), one can characterize systems that has special features 
like inhomogeneities; Nonextensivity is therefore a key role to describe complex systems. 
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